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Abstract 


Second-order optical model solutions to the elastic scattering ampli- 
tude were used to evaluate total , absorption } and abrasion cross sections 
for neutron-nucleus scattering. Improved agreement with experimental 
data for total and absorption cross sections is found when compared with 
first-order (coherent approximation) solutions , especially below several 
hundred MeV. At higher energies , the first- and second-order solutions 
are similar. There are also large differences in abrasion cross-section 
calculations; these differences indicate a crucial role for cluster knock- 
out in the abrasion step. 


I. Introduction 

The assessment of radiation risk to humans, in 
both their natural and space environments, from cos- 
mic radiation is currently an area of active investiga- 
tion. Predictions of biological damage will ultimately 
require a knowledge of the particle fluence spectra 
at the endpoint of interest. In turn, these particle 
fluence spectra are determined from charged-particle 
transport codes, which must contain a description of 
all important physical processes that occur as the in- 
cident ions and subsequent-generation fragment nu- 
clei pass through natural and protective radiation 
shielding. 

Recent studies (refs. 1 and 2) have shown the im- 
portance of using accurate energy-dependent, 
nuclear-interaction cross sections for the determina- 
tion of fluence spectra for high-energy nuclei. A theo- 
retical model for the prediction of cross sections is ex- 
tremely useful, as it cannot be expected that enough 
experiments will be performed for all the collision 
pairs and energies of interest in cosmic-ray studies. 
In previous studies (refs. 3 to 6), a nuclear-interaction 
theory based on an effective coupled-channel solu- 
tion to the nuclear scattering problem was consid- 
ered for the prediction of interaction cross sections. 
A first-order approximation to the elastic amplitude 
was applied within the eikonal approximation for the 
prediction of total, absorption, and abrasion cross 
sections and showed good agreement with available 
experimental data. 

More recently, a second-order solution to the 
eikonal coupled-channels (EGC) model (refs. 7 
to 9) was developed and was found to give im- 
proved accuracy over the first-order solutions in lim- 
ited studies for several collision pairs and energies. In 
this report, neutron-nucleus total, absorption, and 
abrasion cross sections are considered by using the 
second-order ECC model in an effort to improve 
the prediction of nuclear-interaction cross sections at 
cosmic-ray energies. 


In the first-order (coherent approximation) ECC 
model, all nuclear excitations are neglected in the 
evaluation of the elastic amplitude. The second-order 
solution is found in terms of the bordered interac- 
tion matrix, which includes all scatterings from the 
ground to excited states and subsequent decay to 
the ground state. However, cascades between excited 
states are neglected. The second-order solution can 
be expressed in terms of the two-particle form fac- 
tors of the nuclear ground-state, which must include 
the effects of short-range correlations to accurately 
represent the large momentum-transfer region of the 
scattering. In this work, we use a simple model for 
the two-particle form factors based on the nuclear 
matter approximation. Following the procedure in 
reference 10, the combined effects of Pauli correla- 
tions and short-range dynamic correlations are taken 
into account through the use of an effective correla- 
tion length. To simplify the calculations, the zero- 
range approximation is used to evaluate the correla- 
tion terms. The effects of correlations between the 
projectile and target, to first order, are included fol- 
lowing the method developed in reference 6. 

The abrasion cross-section model is based on a 
geometric probability interpretation of the absorp- 
tion cross section. The effects of target correla- 
tions on the abrasion cross section are also consid- 
ered herein by redefining the geometric probability 
for nucleon removal in a manner consistent with the 
second-order ECC solution. Large differences are 
found in this manner between the first- and second- 
order solutions at the lowest energies considered in 
this report. 

Section III contains details of the formalism of the 
second-order solution for the elastic amplitude and 
introduces the second-order model for the abrasion 
cross section. The physical inputs used in our calcu- 
lations are also described. Comparisons are made in 
section IV between the first- and second-order solu- 
tions, and comparisons are made with experimental 


data for calculations of absorption and abrasion cross 
sections. 


II. Symbols 


A f 

fragment mass number 

Ap 

mass number of projectile nucleus 

Ap 

mass number of target nucleus 

B 

slope parameter 

b 

impact parameter 

C{ q) 

correlation function 

fO 

projectile i-particle form factor 

/ 

nucleus- nucleus scattering amplitude 

Inn 

nucleon-nucleon scattering amplitude 


target z-particle form factor 

k 

relative wave number 

k NN 

two-body relative wave number 

k 

correlation length, 0.86 fm 

m 

number of abraded nucleons 

p(b) 

probability of not removing a nucleon 

q 

momentum transfer vector 

r 

internal nuclear coordinate 

T 

kinetic energy of neutron, MeV 

a 

ratio of real part to imaginary part of 
forward two-body amplitude 

r 

second-order eikonal phase 

PCH 

target charge density 

Po 

normalization of density 

<7 

cross section 

X 

first-order eikonal phase 

fi 

solid-angle phase 

Subscripts and superscripts: 

abs 

absorption 

c 

correlation 

dir 

direct 

ex 

exchange 

el 

elastic 


rn abrasion 

NN nucleon-nucleon 

P projectile 

T target 

tot total 

III. Elastic Channel 

In the bordered interaction solution to the optical- 
model coupled-channel equations, the elastic ampli- 
tude is found (refs. 7 to 9) as 

f (q) = I d 2 bexp (— iq ■ b) (exp [»x(b)] cosT (b) - 1} 

( 1 ) 

where k is the relative wave number in the overall 
center- of- mass frame, q is the momentum transfer, 
and b is the impact parameter. In equation (1), \ 
is the first-order eikonal phase, which represents the 
elastic matrix element, and T is the correlation phase, 
which represents the summation over all double scat- 
terings where states are excited from and then de- 
excited to the ground state. The first-order eikonal 
phase is written as 

x(b) = Xdir(b) - Xex(b) (2) 

where the exchange term takes into account corre- 
lation effects between projectile and target nucleons 
(ref. 6). These terms are written as 

Xdir( b ) = J d 2 qexp{iqh)F^\-q)G {l Xq)f NN { q) (3) 

and 

Xex( b ) = 2^J N J d2 9exp(iq b )F^^(-q)G (1 ^(q) 

x J d 2 q exp{iq b)f NN (q^-q f )C{q) (4) 

where F ^ and G ^ are projectile and target ground- 
state one-body form factors, respectively, and 
is the two-body amplitude parameterized as 

/jVA r (q) = - ^NN ex P ^ 

where kjyiy is the relative wave number in the two- 
body system, a is the two-body cross section, B is 
the slope parameter, and a is the ratio of the real 
part to the imaginary part of the forward, two-body 
amplitude. Values for the energy-dependent <r, B, 


2 


and a are found in references 4 and 5. The correlation 
factor is found as 

C(q) = \ G) exp (“ 9 2 / 4 rf 2 ) ( 6 ) 

in reference 6 with d — 1.85 fin -1 . 

The one-body form factor is written in terms of 
the charge form factor as 


where is the ground-state two-body form fac- 
tor of the target. To simplify the evaluation of 
equation (12), the zero-range approximation (refs. 12 
and 13) is used, where 


r 2 (b) = -2 


T 2irA T fNN(0) 
L &NN 


iv: 


,< b > 


(13) 

and where the correlation length l c is 0.86 fm, and 2 
is the z-component of r. 


F^(q) = F CH (q)/F p (q) (7) 

where F p is the proton form factor. For a harmonic 
well distribution, 


The total cross section is found from the elastic 
amplitude by using the optical theorem as follows: 

<Rot = Y lUl ^ = 0 ) ( 14 ) 


Fch = (l ~ s <7 2 ) exp (-a? 2 ) 


(8) Equations (1) and (14) show that 


where values for parameters s and a are from refer- 
ence 5. For nuclei where a Woods-Saxon density is 
appropriate ( A x > 20), 


(r) _ Po 

PCH ^ + exp [(r — R)/c] 


vtot - 4tt J bdb 1 1 - i exp [-Im(x + T)] cos [Re(* + r)] 

- iexp[-Im(x-r)]cos[Re(x-r)]| (15) 

The total absorption cross section is found as follows: 


An exact Fourier transform to obtain the charge form 
factor may be found in a series solution (ref. (11)) 

A'jr 

FcHiq) = — po<fr(q) (io) 


^tot — ^abs F ®el (10) 

where a e i is the total elastic cross section. Integrat- 
ing equation (1) by using dQ « d 2 q/k 2 and by using 
equations (15) and (16) yields 


where 


4>(q) 


= ix Rc 


{ 


— cos (Rq) 
sinh(7 xcq) 


+ 


7 rc sin (Rq) coth(7rcg) 
R sinh(7rc^) 


2c 

7XR 


co 


E(-b 


mcqexp(— mR/c) 1 
(cg) 2 + m 2 ] 2 J 


( 11 ) 


The series in equation (11) converges rapidly, and 
the first three or four terms are accurate for most 
applications. Values for the parameters c and R are 
taken from reference 5. 

The second-order phase T was defined in refer- 
ences 7 to 9 and, for nucleon-nucleus scattering, re- 
duces to 


( 1 

o'abs = 2?r J ^ bdb^l - - exp(-2Imx) 

x |cosh(2Imr) + cos(2ReT) 

All the results of this section (i.e., eqs. (1), (15), and 
(17)) reduce to the coherent approximation in the 
limit T — > 0. 

Another quantity of interest is the abrasion cross 
section, which gives the probability of knockout of 
any number of nucleons. As in references 4, 14, and 
15, the absorption cross section is written as 

a ah , = 2n bdb{l-\P{h)) A ^} (18) 



r 2 (b) = A t ( 2nk NN ) JJ rf2< 7 d V ex p(iq • b) exp(iq • b') 

X /Aw(q)//VJv(q # ) j^-^TG^^(q)G^^(q') 


where P( b) is the probability, as a function of impact 
parameter, of not removing a target nucleon. For the 
second-order solution, it follows from equations (17) 
and (18) that 


+ {A T - l)G (2) (q,q')j 


(12) ^(b)]^ 7, = ^ exp(—2Imx) [cosh(2Imr) -f cos(2ReF)] (19) 
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and the cross section for removal of m target nucleons 
is written as 

am= { A m) 27 T f™ hdb l l ~ ^(b)] m [P(b)]^" m 

(20) 

Because of the dependence of on /nn{ 0)? equa- 
tions (17) and (20) depend on values of a. The 
analogous coherent model solutions are independent 
of a. 

Of special note is the question of whether the 
abrasion model as represented by equation (20) is 
valid when correlations among target nucleons are 
present. The abrasion model is essentially a “clas- 
sical geometric” picture of nucleon knockout. The 
original model is closely related to an independent 
particle assumption (ref. 16) , where cluster knockout 
is not considered. Correlation effects lead to such 
cluster knockouts. (See ref. 17.) Equation (20) does 
not address cluster effects and assumes only that nu- 
cleons are removed. Numerical calculations are con- 
sidered next. Deviations between first- and second- 
order abrasion cross sections suggest regions where 
the abrasion model as represented by equation (20) 
is questionable. 

IV. Results 

Comparisons between calculations and experi- 
mental data are shown in figures 1 to 6 for the total 
and absorption cross section as a function of neu- 
tron energy. The solid line is the second-order model, 
the dashed line is the first-order model, and the ex- 
perimental data (error bars) are from reference 18. 
Above a few hundred MeV, the calculations give sim- 
ilar results; below 300 MeV, the second-order solu- 
tions show better agreement with the data. The rep- 
resentation of T as given in equation (13) may not be 
valid below 100 MeV because of the closure assump- 
tion that is invoked (ref. 8). 

The differences observed at low energies between 
first- and second-order solutions are large for abra- 
sion cross-section calculations. (See figs. 7 to 9.) Dif- 
ferences of a factor of 2 or greater vary with the 
number of nucleons that are abraded. Results for 
abrasion cross sections at 100, 200, and 300 MeV 
are showm in tables 1 to 3. Clearly, redefining P( b) 
to be consistent with the second-order optical model 
leads to a drastically different partitioning of abraded 
nucleons than with the previously studied coherent 
approximation. 

V. Concluding Remarks 

The second-order solution to eikonal coupled- 
channel equations leads to improved predictions of 


total and absorption cross sections for neutron- 
nucleus scattering. Using an abrasion probability 
function consistent with the second-order elastic am- 
plitude leads to dramatically altered predictions of 
abrasion cross sections below r 300 MeV neutron en- 
ergy. The importance of future study of cluster 
knockouts in the abrasion step is clearly demon- 
strated. It would be interesting to extend the cal- 
culations presented herein to other projectiles. Cal- 
culations with ablation effects should be made and 
compared with experimental data for fragmentation 
cross sections. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November 15, 1991 
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Table 2. Abrasion Cross Sections for Neutron 27 A1 
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(c) Four-nucleon removal. 
Figure 7. Concluded. 
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(c) Three-nucleon removal. 
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(d) Four-nucleon removal. 
Figure 9. Concluded. 
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